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We calculate the bulk contribution for the doubly degenerated largest eigenvalue of the transfer 
matrix of the eight vertex model with an odd number of lattice sites N in the disordered regime 
using the generic equation for roots proposed by Fabricius and McCoy in We show as expected 
that in the thermodynamic limit the result coincides with the one in the N even case. 



I. INTRODUCTION 



In the last years it has been seen that in exactly solvable systems there are interesting features which distinguish 
between systems with an an odd and even number of lattice sites 1 II S, 0, 0, E 0, Ed, Ell ■ The properties of finite 
size systems are qualitatively different in the two cases and also totally new features occur in the odd case. 

Recently in a series of papers [H,[lH Fabricius and McCoy have investigated the properties of the Q72(v) matrix 
introduced by Baxter in Appendix C of [l3| in his study of the eight vertex model (for more information on the eight 
vertex model see [14j). This matrix is defined for systems at the root of unity 



mK 

v = — (1) 

where K is the complete elliptic integral of the first kind of modulus k, m and L are integers and the zeroes of the 
eigenvalues of Q72(v) play an important role in the calculation of the eigenvalues of the transfer matrix of the eight 
vertex model. The result of this investigation showed that the distribution of these zeroes in the case of an even 
number of lattice sites N is qualitatively different from the N odd case. In the N even case these zeroes appear either 
as L-strings of the form + 2jK/L,j = 0, • •• ,L — 1 or pairs and + iK' called "Bethe's roots" which satisfy 
"Bethe's equation" 

\h(uf-v)J ~ l} 1 Kuf- U B-2r 1 ) U 

where 115 is the number of pairs, H{v) and 0(f) are elliptic theta functions of modulus k (for definitions and some 
properties see Appendix [A"|) and h(v) = H(v)Q(v). 

In the N odd case when m is odd (this is the case where the matrix Q72(u) exists [l|, [HI) we do not have Bethe 
roots or the L-strings. We do have however the following type of pairing: for every set of roots Uj we have the second 
pair of roots Uj + iK' which also gives an eigenvalue of Q72(v) and these two eigenvalues satisfy the "TQ" equation 
with the same eigenvalue of the transfer matrix T(v). These roots satisfy (3.3) of [l[ 

h(uj -rj)J H{ Uj - u k - 2??) 

This equation differs from "Bethe's equation" in that the right hand side contains the function H(v) and not 
h(v) and the number of terms in the product is N. For the m even case numerical investigations showed that the 
ground state has the same structure of roots and they obey also ([3]). The equation ([3]) is new in the literature and 
the purpose of this paper is to calculate the largest eigenvalue for the transfer matrix with N odd (m even or odd) 
when the largest eigenvalue is doubly degenerate. 
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We make this observation more precise by focusing on the case with N odd which is the main interest of this paper. 
Baxter's [3], [IB] original computation of the eigenvalues of the transfer matrix T(v) of the eight vertex employs the 
so called "TQ" equation 

T(v)Q 72 (v) = [pQ(0)h(v + v)] N Q72(v - 2 V ) + [pQ(0)h(v - rj)] N Qn(y + H (4) 

where Q72 (v) is an auxiliary matrix with the properties 

[T(v), Q 72 (v')] = [Q72(w), Q 72 («')] = (5) 

and 77 is restricted at the root of unity condition ([I]) . Also in [l3| it has been showed that if we parametrize the vertex 
weights in terms of p, v, ?y, k as 

a(v) = pQ(2r})G(v-r))H(v + r}) (6) 



b(v) = pQ(2rj)H(v- rj)e(v + 77) (7) 

c(v) = pH(2r))e(v - r))e(v + 77) (8) 

d(v)=pH(2r))H(v-r])H(v + r)) (9) 
then transfer matrices with different values of the parameter v commute 

[T(v),T(v')]=0 (10) 

The operator S defined as 

N 

s=n^ (11) 

commutes with T(v) and Q.72(v) 

[S,T(t>)] = [S,Q 7a (t;)]=0 (12) 

has eigenvalues (— 1)" = ±1 and under spin inversion a z — > — a z we have S — » (— 1) S . Under spin inversion the 
eigenvalues of the transfer matrix are invariant so in the case of TV odd the eigenvalues of the transfer matrix are 
doubly degenerate one with v' = and one with v' = 1. In the N odd case the matrix Q72(w) which satisfies (1.5) 
and (1.6) of [l[ and the "TQ" equation exists only in the case when 77 is of the form (TT|) with m odd and from [l| we 
know that the eigenvalues of the Q72 matrix can be written in the factorized form 



N 

! 72 {v)=K,{q,u k )e- i ^JlH{v-u k ) (13) 



k=l 

with JC{q, Uk) a constant and u k the zeroes of the eigenvalues of Q72- As we have said in this case with m odd and L 
even or odd the roots satisfy (J3j) and we do not have Bethe roots or L-strings. 

In the case of m even with L odd numerical investigations presented in [l| revealed the following new features as 
r) — > mK/L: some of the roots moved to form L-strings giving a contribution which cancels from the generic equation 
of roots ((3|) and the rest arranged in pairs u k and u k + iK' . In this cases (j3]) reduced to the "Bethe's equation" ([2]). 
However the roots of the solutions without the L-strings or paired roots still satisfy ([3|) . Also it should be mentioned 
that in both cases (m even or odd) the ground state is characterized by solutions with N roots satisfying ([3]). 

These roots of the eigenvalues of the Q72(f) are the main ingredient in calculating the eigenvalues of the transfer 
matrix. The free energy per site of the infinite lattice is related to the largest eigenvalue of the transfer matrix by 

/ = -k B T lim iV- 1 In \T max (N)\ (14) 

where ks is the Boltzmann constant. 

Originally in 1971 Baxter [l5[ calculated the largest eigenvalues for the N even case using a new perturbation- 
expansion technique. In 1973 Johnson, Krinski and McCoy [l7| reproduced the result of Baxter and calculated also 
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the next largest eigenvalues using the integral-equation method which was applied previously to other models solvable 
by Bethe ansatz. In 1988 Kliimper and Zittartz [l8[ using an inversion relation for the eigenvalues of the transfer 
matrix obtained all the eigenvalues and in 1989 they classified the excited states f]~Qj j . 

In our paper for the calculation of the largest eigenvalue of the transfer matrix we are going to use the integral 
equation method (Sect. IV of [ljj and Sect. VI of [2Cj|]). The plan of the paper is as follows: in Scction[TT]wc are going 
to obtain and solve the linear integral equation for the density of roots. In Section HTT1 we are going to use the explicit 
expression for the density of roots to calculate the largest eigenvalue of the transfer matrix in the thermodynamic 
limit which is the main result of this paper. Using the result of Baxter obtained for the N even case (for a 
sketch of the derivation in the disordered region see [2l|) we will show that the bulk contributions are the same. The 
technical details of the calculations are relegated to the Appendixes. Our analysis will be restricted in the disordered 
regime; more precisely in the region 

< q = e- 2 ^ < 1 (15) 

< 7} < K (16) 

with K > and K' > where K is the complete elliptic integral of the first kind of modulus k and K' the same 
integral of complementary modulus. Also we are employing the convention that the theta- functions whose modulus 
is not explicit should be understood as k. 

II. INTEGRAL EQUATION FOR THE DENSITY OF ROOTS FOR N ODD 

The main feature of the eight vertex model with N odd is the fact that the largest eigenvalue of the transfer matrix 
is doubly degenerate. In order to calculate it (for both m even and odd) we are going to use © the generic equation 
of roots proposed by Fabricius and McCoy in [l[ and the following set of roots (see Table I in |l|) 

u) = ivj + K with - K' < v j < K' (17) 
iK' (18) 



.2 _ „,1 



Using the sum rules obtained in [1] 



K 

3=1 



even integer (19) 



A' 



1 777 [U i 

~= even integer — v 1 — JV (20) 



K 

3=1 



we have for the first set of roots {u]}, v = 0, v' = 1 and for the second set of roots v = N, v' = 0. 

Using (|A7[) and (|A5|) eq. [3] for the largest eigenvalue can be written for both sets of roots as 

Hi(ivj + rj)&i(ivj +rj) \ N yj H(ivj — ivk + 2r/) 



H\{—ivj + r])®i(—ivj + rf) J ±± H(—ivj + iv k + 2rj) 



TT ^3-^-r^U (J = l r ..,N) (21) 



Define the functions 



Hi(p$ + ix)Qx(p£ + ix) 



Mx,v=^ 7r2 , zcr, , ( 22 ) 



(23) 
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Both of these functions are analytic in a strip around the real axis (see Appendix [B| defined by the conditions 

\Tmx\ + \p£\<K (24) 

and 

\Xmx\ + \(p£-K)\ <K (25) 
respectively. In terms of these functions (|2ip takes the form 

N 

NA 1 (v J ,i 1 ) = 2TrI J +Y / B 2 (v 3 -v k ,7 1 ) (j = l,...,N) (26) 

k=l 

where Ij are integers specifying the branches of the logarithm. For the largest eigenvalue we have Ij — Ij—i = — 1, J — 
2, • • • , TV. As TV — > oo the number of roots becomes infinite and their imaginary parts fill the interval [-K 1 , K']. The 
number of roots with the imaginary parts in the interval [x, x + dx] is given by Npo(x)dx where po(x) is the density 
of roots. Using 

1 - r K ' 

jT7 2 /(«*)=/ f(x)po(x)dx, N ^ 00 (27) 

fe=l J-K' 

in the large TV limit the coupled nonlinear equations are replaced by the linear integral equation 

A[{x,j]) = -2ix Po (x) + B' 2 (x- x') Po {x')dx' (28) 



where the prime denotes the derivative with respect to x. This equation can be solved using the Fourier transform 
(see Appendix ICj) obtaining 

^ oo 

Po(x)=w J2 Po(n)e-^ (29) 



2K> 

71— — OO 



with 



Pb(0) - 1 

po(n = odd) = (30) 

i 

po{n = even) 



ooBh[n(^-S0] 



III. LARGEST EIGENVALUE OF THE TRANSFER MATRIX 

It can be shown (see Appendix [Pi for a proof in the domain v + 2r\ < K , < v, rf) that in the "TQ" equation 

Tno{v,v) = [pQ(0)h{v + r?) — — ; — + ^6(0)/i(u - 77) ^ — 

= T+ (v,v') +T NO {v.v') (31) 

T^ (v, v') exponentially dominates T^ (v, v') for large TV so we will focus our attention on the first term. From [l[ 
we know that the eigenvalues of the Q72 matrix can be written in the factorized form 

TV 

Q 72 (v, v') = K{q, u k )e-^ Yl H(v - u k ) (32) 

fc=i 

with K,(q, Uk) a constant and u k the zeroes of the eigenvalues of Q72. Explicitly we have for T^ (v, v') the following 
expressions for the first set of roots (fl"7|) when v = 0, v' = 1 

N 

TUv, v' = 1) = [ P mh(v + V)] N II ^V + g^ (33) 

^1 H(iv k + K - v) 
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and for the second set of roots (fT8|) when v = N, v' = 0, using (|A7 



r+ o(Ul ✓ = o) = [pemiv + ,r n (34) 

fe=i Q(wfc + -fv - v) 
After we have used (|A5|) and the Jacobi transformation (|B2[) and (|B4[) we obtain 

^o"'" )= z 30 ^ + M e k k' — _ — — (35) 

^ £A B(wfe + iv + (1 — v')K, k') 

Defining 



3o (x,,)=ln(e k k' Q{x + w + {1 _ vl)K ^ kl) ) (36) 



and taking the logarithm of (|35[) we have 



iV 



lnT+ ( W , t/) = N\n[ P e(0)h(v + 9o(^k, v 1 ) (37) 



fe=i 

As in the previous section the sum in (f3"T)) can be transformed in 



k=i 

or in terms of Fourier coefficients 



Y>S(u fe y)=iv/ dis5(sy)po(a!) (38) 
J—K' 



N N +oo 

Xl^oK^') = ^7 XI 9o( n > v ')Po(-n>) (39) 

fc=l n= — oo 

It follows from (jAllj) and (|A12I) that 

,. irf-nv — n 2 — iiix) 2(— l)[( 1 ~ l/ )™1 . { nitix + iv — in) \ . / nirir]\ , , . 

g±(x, v 1 ) = — ! — + > — — ^ — ^ sin — sin (40) 

y ° y ' ; KK' ^nsinh(^) V K' ) \ K' J v ' 

valid for 

\lm(x + iv - irj) \ + r) < K (41) 

so the Fourier coefficients are 



2if' / »,(.-,) (-DK^^^sinhf^P) n 

I ^ rw x — 1_ 



g£(n,v') = — I e^^ K —± S ' - -^cos(wr) I (42) 

• 1 sinh(iff) K V '< 



n 



Collecting all the results and noting the fact the Fourier coefficients of the density of roots are zero for n odd we 
obtain for the largest eigenvalue of the transfer matrix (for both sets of roots) 



g sinh[2n( F -ff)]sinh(^ 

n=l 



ns inh(22H^)cosh [2n(^ - §£)] 



(43) 
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in the domain 

< r\ < v < K (44) 



Even if our result was obtained in the domain (|44[) we see that we can consider the analytical continuation of (I43|) 
which is valid in 

<7] < K, i] < v < 2K - 77 (45) 

Now we will show that our result for the largest eig envalue of the transfer matrix in the disorder regime (|43|) coincides 
with the result obtained for even N (A. 51 of [21j). The N even result obtained for the normalization factor p' — 
2iJ 1 (0| g 1 / 2 )- 1 9(0| 9 )- 1 is 



1 2»Vr«' r/1 00 sinhfg^sinhf 2 "^-^ ) 
1 lnT+ E (v') = hyOriv' + V '\q') + ^ {V V ] + V ) f (46) 



" = 2* W= 2T T= 2if 9=6 H{x)=e l {—\ q ) (47) 



AT 

with 



It can be easily seen that (|46[) and (j43|) are equal if we consider in our result the same normalization factor and the 
identity (jA"T0l) . 

IV. DISCUSSION AND CONCLUSIONS 

We have calculated the doubly degenerated largest eigenvalue of the transfer matrix of the eight vertex model with 
an odd number of lattice sites in the thermodynamic limit. The next interesting problem would be the calculation of 
the low lying excited states for the odd case in a manner similar with the computations in (l7| . In the m odd case the 
calculations should be accessible, however in the m even and L odd case an interesting problem occurs. As we have 
said in the m even and L odd case some of the roots move to form L-strings which drop from the generic equation 
for roots ([3]) while the rest arrange themselves in pairs. For these states the usual procedure of calculating the 
energy of the excitations using linear integral equations cannot be applied. This means that in this case the excitation 
spectrum has two different kind of states. The m even and L odd case is different because in this case some of the 
eigenvalues develop extra degeneracies beyond the double degeneracy present in the odd case. The 77 = 2K/3 case is 
particularly interesting due to the fact that in this case it is believed that the degenerated largest eigenvalue is exact 
even for a finite lattice 

^ In T NO (v)= a + 6, 7? = 2K/3 (48) 

where a and b are the Boltzman weights of the model (see formulae |B] and [7]) . This was advocated in (23l . [24], [25| and 
verified analytically for TV up to 15 in [2lj . 
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APPENDIX A: ELLIPTIC THETA FUNCTIONS 

Let K be the complete elliptic integral of the first kind of modulus k , K' the same integral of the complementary 
modulus k' = (1 — A; 2 ) 1 / 2 and q — e~* K ' K . The Jacobi theta-functions of nome q are defined by the relations 

00 

H(x) = 2gV4 sin J[ (l - 2 g 2 » cos (^) + g 4 ") (1 - , 2 ») (Al) 

n=l 
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valid for \Imx\ + \Im(y - iK')\ < K' 



Hjjx + y) = ( H{x + y-K 



ittx \ 2(— l) n . fnitx\ . f mr(y — iK') \ 

> , sin ( — — sin — ^— '- (A14) 



H^x-y) V H(x-y + K) 

- 2(- 

K ' ^ nsinh(s^l) V X 7 V X 

n— 1 \ K / v 

valid for |2mx| + |Zm(y - iK')\ < K' 

APPENDIX B: ANALYSIS OF THE RELEVANT FUNCTIONS 

In order to facilitate our analysis of A p {x,^) and B p (x,^) functions we are going to perform Jacobi transforms on 
the theta-functions involved. The relevant formulae are (Chap. XV of [14] or Chap. XXI of (H) 



H(ix, k) = iy/K' 'jKeSwH{x, k') (Bl) 



Hi(ix,k) = y/K'/Ke^Q(x,k') (B2) 



S(ix, k) = \j K 1 1 Ke^KK 7 ' Hi (x, k') (B3) 



Qi(ix,k) = y/K'/KeZKfcQifak') (B4) 

where we have explicitly expressed the modulus dependence of the theta-functions involved. 
After the Jacobi transformation and using (|All|) and (|A12j) we obtain for A p 



Ap(x,£) = ihx 



Hi (p£ — ix, k)&i(p£ — ix, k) 
Hi {pS, + ix, k)Oi(p£ + ix, k) 



i 2»^ P ;x Qiip^ + x, k')Qi(ip£ + x, k') 
i m e kk> 



®(-ip£ + x, k')Qi(-ip£ + x, k') 

2npfx 2i ( 2mrx\ ( 2niripf\ ,„ , 

= — + > To — ttt sin sin — (B5) 

KK> ^ nsm h(%^) \ K> ) \ K> ) y ' 

where the expansion is valid for \Xmx\ + \Tm(ip£)\ < K and we specify the branches of the logarithms by also 
demanding that no cuts of A p cross the real x axis and A p (0, £) = 0. 

Analogously after the Jacobi transformation and using (|A13[) we obtain for B p 



( H(p£-ix,k) 
\H(p^ + ix,k) / 

vi* H{ip^ + x, k') 



= iln (H^ H (-ipt + x,k>)) 



, P€\ kx 2i . / tlttx \ . fnw{ip^-iK)\ 
1 - k ) K> + g nsinh(^) Sm ^ { K> J (B6) 



where the expansion is valid for \Tmx\ + \Tm(ip£, — iK)\ < K and we specify the branches of the logarithms by 
demanding that no cuts of B p cross the real x axis and B p (0, £) = 0. 
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APPENDIX C: DENSITY OF ROOTS 



We define the Fourier coefficients of a function periodic (quasiperiodic) of period (quasiperiod) 2K' by 

r K ' 

f(n) = / dxf{x)e-i^ (CI) 



>-K' 

In terms of this coefficients the original function is expressed as 



oo 



2if 

n— — cx) 



After the Fourier transform the linear integral equation for the roots takes the form 

A[(n,ri) = -2wpo(n) + B' 2 (n,r])po(n) (C3) 



po{x) = J- y _%^L e -^ 



(C4) 



In order to obtain explicit expressions for the density of roots we need to calculate the Fourier coefficients of B' 2 
and A[ . From (|B6[) we obtain the following expansion 

/ 2ri\ irx 2i /mrx\ ( wr(2wj — iK) \ ,„ , 

*(*.*> - (l - i) F + £ — 5^ ™ (-gr) sin ( ^ ' ) (C5) 

so the Fourier coefficients for £?2 are 

^(n,„) = 2 / mM ri?: *' jJ (C6) 



sinh 1 



In the case of A[ using (|B5[) we have the following expansion 



. . , 2irrix 2i . / 2n7rx\ . / 2n-n:iri\ 

Ai fx, 77) = V- + > 7^ — 77T- sm sin (C7) 

V " KK< ^nsinh(22f^) \ K> ) \ K' ) K 1 



so the Fourier coefficients are 

AM = 



l^(n = odd, 77) = (C8) 

■2>/ s sinh(^) 

A\ (n = even, 77) = -47r ,v 



sinh(^) 
Using ([C6]) and JCS]) and the identity 

sinh a + sinh(26 — a) = 2 sinh 6 cosh(6 — a) (C9) 

we obtain 

Pb(0) = 1 

po(n = odd) = (CIO) 



APPENDIX D: RELEVANT TERM IN THE "TQ" EQUATION 

The logarithm of the ratio of the two terms in the "TQ" equation is 

ln ihM = Nlxi ( h(y + v) \ ln ( Qnoiv-tn) 



T NO (v) \h(v- V )J \Q 7 2o(v + 2rj) , 

Let's analyze the first term in the RHS of (|D1[) . After the Jacobi transformations it becomes 



Nln 

so using (|A13|) and (|A14[) we obtain 



h(v + ri)\ _ ri / H(iv + irj, fc')ifi (iv + iri. fc') 

I = iv in e kk' 

h(v — rj) J \ H(iv — irj, k')H\{iv — irj, k') 



/h(v + n)\ 2nvn ( 2niv\ ^ 2N ( '2mrn\ f 2rm(K 
Win 77 -r = -N——!- + iN -2tt - —— + V .„ y . sinh — — * sinh ' V 



h(v-r))J KK' V K J ^nsinh(2^«) V X' 



where the sum is a positive number (0 < r\ < v < K). 

For the second term in the RHS of (|D1[) using (f3"2")) . (|A5|) and the Jacobi transformation we obtain 

lu ( Qnojv -2 V ) \ = ^ / H 1 (iv j -v + 2r}) 
11 \Q72o{v + 2 V ) J ^ i n \H 1 (iv j -v-2r ] ) 

N 



Defining 



the sum can be written as 



E, / i-xivivj+iv) Q( Vj + iv - 2tn 7 k') 
ln e kk' — 
\ Q(vj + iv + 2it],k') 



. 2ni V ( : c+iv) Q(x + iv — 2in, k') 

q(x) = In e w 1 — 

' 8 (a; + iv + 2iri, k') 



N I + oc 

fc=l \ n= — cx 



5(n)p (-n) + 0(l) 



From (|A11|) we have the following expansion 



. . 27TW77 2iirnx v 2 / rwrfa: + it?) \ / 2irnrn\ 

g(x) = 7 + > -. — rrr- sin — - sin 

yv ; KK' ^nsinh(2ff) V / \ K' ) 



so we obtain 



AT 



V \-N\ ^ sinh (jga) 1 

tr U^ + r _^ o 2nsinh(^)cosh[2n(^-^)] + "V 



From (1031) and (00911 we see that 



lnSoM = 7Vc+C(i) 
TnoW N 



where c is a number whose real part is positive. 
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